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Abstract Clés is an évolution code recently developed to 
produce stellar models meeting the spécifie requirements of 
studies in asteroseismology. It offers the users a lot of choices 
in the input physics they want in their models and its versa- 
tility allows them to tailor the code to their needs and imple- 
ment easily new features. We describe the features imple- 
mented in the current version of the code and the techniques 
used to solve the équations of stellar structure and évolution. 
A brief account is given of the use of the program and of a 
solar calibration realized with it. 

Keywords stars • stellar structure • stellar évolution 
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1 Introduction 

Since the early 70s, the asteroseismology group at the In- 
stitute of Astrophysics and Geophysics of Liège has been 
using an évolution code derived from the Henyey code. This 
code has been continuously updated. Nevertheless, with the 
progress of asteroseismology, it became clear that the fre- 
quencies and the stability of the oscillation modes were ex- 
tremely sensitive to détails of the model, which were unim- 
portant for the computation of stellar évolution. A few years 
ago, it was thus decided to write a new code, meeting the 
spécifie requirements of our studies in asteroseismology. This 
code has been named Clés, the acronym of Code Liégeois 
d'Evolution Stellaire. It is still in an active phase of develop- 
ment. A number of persons have contributed and contribute 
to its development. The main objective of this effort is to 
have a versatile stellar évolution code producing models pré- 
cise enough to be useful for studies in asteroseismology. We 
try to maintain a clear structure so that every user will be 
able to tailor the code to his own needs. 
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Clés has been developed on a GNU/Linux system (more 
precisely a Slackware one). No effort has been made to port 
the software to another platform. It is written in FORTRAN 
77. A few relatively standard extensions of the language pro- 
vided by the GNU Fortran compiler g77 have been used. A 
few scripts use the Bash shell and the script building the 
Makefile is a Tel script. 

The présent version of the code is numbered 18. The 
code is not unique. As it is easy to customize, several ver- 
sions of the code have been developed from the standard ver- 
sion, implementing différent features needed by the works 
already realized or in progress. The most interesting of thèse 
will be included in the next standard version of the code so 
as to préserve its uniqueness. The features described below 
do not necessary belong to the standard version. After a brief 
reminder of the équations we have to solve, we describe the 
physical assumptions. Two sections are then devoted to tech- 
nical aspects and we end with considérations on the use of 
Clés, a solar calibration and a brief discussion. 

Due to limitations in the physics already implemented, 
the présent version of Clés is not intended for the hélium 
burning phases of the évolution. 



2 Structure and évolution équations 

2.1 Structure 

The quasi-equilibrium structure of a star is governed by the 
following équations (we use well-known notations): 
dp Gmp 



— =4nr^p, 
dr 

— =Anr p (e- 
dr 

dT_ 
dr 



GinpT 
~r^ 



(1) 



(2) 



(3) 



(4) 



In Q, the luminosity L does not include the neutrino lu- 
minosity and the term e dénotes the rate of production of 
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nuclear energy from which we have already excluded the 
neutrinos losses. In the same équation, the term deserves 
a comment. It hides a time derivative and governs the évo- 
lution rate when no significant nuclear energy sources are 
available, for instance in the beginning of the pre-main sé- 
quence phase. It reads 



du P dp dS ^ dm 

= = ~T — y Hi — -. 

dt p2 rit 



dt 



dt 



dt 



(5) 



where /i, is the chemical potential of élément i and n, the 
number of moles per gram. In t his expression, some peo- 
ple keep only the term —TdS/dt dKippenhahn and Weigerj 
[l990), but this neglects the effects of changes in chemical 
composition. An enlig htening discussion of the correct form 
of this term is given by IStrittmatter et afl (Il970l) . 

In S, V stands for d\nP/d\aT and is computed dif- 
ferently according to whether the considered point belongs 
to a radiative or convective zone. In a radiative zone, the 
diffusion approximation allows us to write the température 
gradient as 



V = Vrad = 



3kLP 



IGnacGinT'^ 



(6) 



The expression of V in a convective zone is discussed in the 
section devoted to convection. 



do not dépend on r and (|7]| is replaced by its integrated form 
over the entire convective zone: 
dXi r'«2 



dt 



m] 



Y^Rijdm 



(8) 



where nie is the mass of the convective zone and indexes 1 
and 2 refer to the bottom and the top of the convective zone 
respectively. 

Initial conditions. The first model of a séquence is always 
on the Hayashi track. It is chemically homogeneous and its 
central température is between 4x10^ and 5 x 10^ K. The 
initial chemical composition is defined by X and Z, the mass 
fractions of hydrogen and metals. The métal mixture and the 
isotopic abundances can be changed by the user. It is how- 
ever his responsibility to select appropriate opacity tables. 

The standar d version of Clés comes w ith two mét al mix- 

tures, GN93 dGrevesse and Noels ^ [I993t and AGS05 dAsplund et al.[ 
l2005h and the ' Anders and Grevasse! isotopic abun- 

dances. It must be stressed that for Li, Be and B, contrary 
to other éléments, we do not use solar photospheric abun- 
dances. As thèse éléments may have undergone nuclear pro- 
cessing at the basis of the solar convective zone, we have 
adopted their meteoritic abundances. 



Boundary conditions. The differential équations of structure 
must be supplemented with boundary conditions. The condi- 
tions at the center are obvions, m and L must vanish as P. In 
the external layers of the star, the diffusion approximation 
progressively breaks down as the optical depth approaches 
unity. Rather than integrating the full radiative problem, it 
is usual to join smoothly the interior solution to a precom- 
puted atmosphère model, which imposes, in this way, exter- 
nal boundary conditions on the interior model. 



2.2 Chemical évolution 

The équations governing the rate of change of the abun- 
dances of the éléments may be written in the form 



dXi 
dt 



J_d_ 



pXjWi 



(7) 



Rij expresses the contribution of the j-th nuclear reaction to 
the variation of the species i. It dépends on the abundances 
of several species in a non linear way and needs no further 
explanation. w, is the diffusion velocity of élément /. The 
computation of this term is described in section [331 

In convective and overshooting zones, the différent 
species are rapidly mixed. In a few évolution codes this rapid 
mixing is described as a diffusive process with appropriate 
coefficients. In Clés, convective mixing is an instantaneous 
process and the chemical homogeneity of convective zones 
is imposed. In a given convective zone, the abundances X, 



3 Input physics 

3.1 Equation of state 

We do not compute the équation of state (EOS) in the évo- 
lution code. A few EOS are provided in the form of pre- 
computed tables which are interpolated. We are particularly 
careful with respect to the consistency between the différ- 
ent thermodynamic quantities provided by the équation of 
state. The independent variables are logp, logT, X and Z 
(log stands for the décimal logarithm). At each grid point, 
we store only logP, logC,., Pp = {dlogP/dlogp)T and Pt = 
{d logP/d log T)p for the gas. In the EOS routine, thèse quan- 
tities are interpolated, the easily computed contribution of 
the radiation is taken into account and ail other thermody- 
namic quantities are deduced from well-known identities. 
The contributions of the radiation are included as follows. 



P = Pgas + Prad with Prad = \aT'^ . 



AaT^ 



gas 



Pp=PPp.gas with l5=Pgas/P, 
PT=pPT.gas+ 4(1-15). 

Adiabatic exponents and Cp are then computed, 
PtP 



r3-i = 



C,pT ' 



ri=Pp + (r3-i)PT., 



(9) 
(10) 

(11) 

(12) 



(13) 



(14) 
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Cp = FiCv/Pp ■ 



(15) 



Three EPS tabl es have been use d with Clés impl ement- 
ing CEFF dChriste nsen-Dals gaard a nd Dapped . 1992h . OPAL 
2001 and OPAL 2005 (Rog ers and Nayfonovl . l2002h . We have 
computed the CEFF table with a routine kindly provided by 
J. Christensen-Dalsgaard. A métal mixture defined by the 
abundances of ten metals can be input to the routine. As 
the EOS is not very sensitive to the métal mixture, we have 
computed only one table with a mixture as close as possible 
to the GN93 mixture. The OPAL 2001 and 2005 EOS are 
available in tabular form for a fixed métal mixture of four 
éléments (C, N, O and Ne). As the C, from OPAL 2001 ta- 
bles are known to be inaccurate (Boothrovd and Sackmann, 
I2003h . Cy has been computed from the other tabulated quan- 
tifies. 

For the necessity of diffusion a routine computes the 
stage of ionization of ail the éléments involved in the dif- 
fusion process (Saha's équation). 



3.2 Opacity 



The s tandard version of Clés uses OP AL opacities ([I glesias and 
1 19961) . completed with the opacities of lAlexander an d Ferguson 
(Il994) at low température. Both tables have been smoothly 
merged into a single table. In the température domain log T G 
[3.9,4.15] where the opacity is defined in both tables, we use 
an opacity K defined as the average 



log ïf = ( 1 - ) log JCa/t + log KoPAL , 



(16) 



In this expression, 6 is the unique third degree polynomial in 
log T, with vanishing derivatives at both ends of the interval 
and taking zéro and unity values at the lower and higher ends 
respectively. 

A few opacity tables have been built for différent métal 
mixtures. Appropriate tools allow the user to build new ta- 
bles in the format required by Clés from tables generated 
at the OPAL site and the tables of Alexande r and Ferguson. 
For pa rticular studies, tables using the new 'Ferguson et al.l 
(|2005|) opacities and OP opacities (Badnell et al., 2005) have 
also been built ( Miglio et al.LI2007bO . 

For the computation of models with local variations of 
the métal mixture, several tables with différent métal mix- 
tures are loaded and interpolated (M iglio ctïà 

Up to now the effects of conduction on opacities have 
been neglected. 
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CNO cycles: 
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(35) 



We follow thoroughly the combustion of and ^Li. 
Only unstable species (^Be, '^N, '^O, '^F and '*F) are sup- 
posed to be at equilibrium. Though the main reactions of the 
hélium burning phase have already been implemented, we 
have yet to improve our code (semi-convection, équation of 
State, opacity) to be able to accurately follow this phase of 
the évolution. 

The reaction ra tes are computed using the a nalytical ex- 
pressions given by ICaughlan and Fowled (Il988h . For the re 



action ''^N(p,7)^^Ô, we use the cross-section given by lFormicola et al.l 
(120041) . A variant of the program using th e NACRE reaction 
rates in their approximate analytical form dAngulo et ami999h 

has been written to facilitate the comparisons with CESAM 

('Mont albân and LebretôriLl2005l : lMontalbânLl2ÔÔ5l : lMontalbân et al.l 
,2005) . 

Th e screening factors are computed according to lSalpetëJ 
(119541) . 



3.3 Nuclear reactions 

The following reactions are included in our nuclear reaction 
network. 

p-p chains: 

2'H^2H + e+ + v (17) 
^H+^H^^He+y (18) 
2^He^'*He + 2 (19) 



3.4 Atmosphère 

At the last point of a stellar model Clés realizes a smooth 
junction with a pre-computed model atmosphère. To that 
end, tables giving the outer boundary conditions in terms 
of g, Tef f , X and Z have been computed from Kurucz atmo- 
sphères dKuruczL Il998h at différent optical depths (photo- 
sphère and T = 1, 10 and 100). We have also prepared tables 
of boundary conditions from radiative gray Eddington atmo- 
sphères. 
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3.5 Convection 

A région of a star is dynamically stable against convective 
movements, if Ledoux's criterion is satisfied: 

dlnp 1 dlnP 



dr 



n dr 



<0. 



(36) 



In a chemically homogeneous zone, it is équivalent to the 
criterion of Schwarzschild: 



ad, (37) 

where Vad = (^3 ^ 1 ) /^l is the adiabatic gradient. In prés- 
ence of a chemical inhomogeneity, both criteria cease to be 
équivalent. In term of V, Ledoux criterion reads 

V<Virf,, (38) 

where V ux dérives from the rewriting of équation ( 1361 ). It is 
sufficient for our purpose to give its expression for a com- 
pletely ionized perfect gas 

Vu,=V„rf + V^, (39) 

where = d\npL/d\nP is the gradient of the mean atomic 
weight. For a non homogeneous région to be in radiative sta- 
bility, it is not enough to satisfy Ledoux's criterion. In that 
case, Schwarzschild's criterion is a necessary condition for 
vibrational stability. In a star, due to nuclear processing, the 
mean atomic weight generally decreases from the centre to 
the surface and the term is positive. We have then 

yad<yUx, (40) 

and when Schwarzschild's criterion is satisfied, Ledoux's 
one is satisfied too. However, when diffusion is taken into 
account, gradients may develop in small régions of the star, 
where this inequality is reversed and Ledoux' criterion must 
be explicitly taken into account. 

We have imp l ement ed the usual mixing-length theory 
of |B o hm- Vi tensd d 19581) . also exposed in the textbooks of 
IÇoxand Giulil dl968h"ândrKippenha hn and Weig ër3 dl990l) . 
It is a local theory, in which the température gradient is ob- 
tained very simply by solving the cubic équation 



-r'+r=A{Wrad-y, 



ad) 



(41) 



where A is expressed in terms of the mixing-length ^, gener- 
ally defined as a multiple of the scale-height, l = aHp, with 
coefficient a chosen by the user. 



PtP 



2PpP 



KCpP' 



\2acT^ 

The real température gradient is then expressed as 

(r + i)v„d 



v = 



r^^ad- 



r + i 



(42) 



(43) 



The variant due to iHenvev et al.1 dl965l) has also been 
used for the description of convection at low optical depth. 

The Full Spectrum Turbulence (FST ) theory of convec- 
tion dCanuto and Mazzitelli, 1991 , 19921) has also been im- 
plemented dMiglio and MontalbânLl2005h . 



No effort has been made to put a grid point at the bound- 
ary of a convective zone. Though the mixing of the material 
inside convective zones takes into account the précise loca- 
tion of their boundaries, this results in a small diffusion of 
numerical origin at the boundaries. 



Overshooting. Overshooting (or undershooting) displaces the 
boundary of a mixed région from to r^y given by 



Tov = rc± Uov mm{Hp,h) , 



(44) 



where h is the thickness of the convective zone and the over- 
shooting coefficient tto,, is chosen by the user. 

In the overshooting région, the température gradient can 
be chosen adiabatic (only for the core), radiative or com- 
pute d according to the prescription of the mixing length the- 
ory dGodartLl2007h . 



Semi-convection. No semi-convection treatment has been im- 
plemented. For low-mass main-sequence stars, this results in 
an uncertainty on the position of the boundaries of the con- 
vective régions, at the limit of the core when it is increas- 
ing and also at the botto m of the convective envelope when 
diffusion is included (Miglio et al.L l2005l : iMontalbân et al.l 
120071) . 



3.6 Diffusion 



We follow the theory of stellar diffusion developed bv iThoul et alj 
(1994). From their fundamental équations (12), (13), (18) 
and (19), we eliminate the électron density and express the 
diffusion velocity w, in terms of the gradients of abundances 
of ail ions, of pressure and température. 



V. dXj 



-ap- 



dlnP 



dr 



d\nT 



dr 



(45) 



Though it is possible to eliminate one of the ion abundances, 
it seems that keeping ail the abundances ensures a better sta- 
bility to the solving process. 

In the standard version of the code, three groups of él- 
éments are distinguished for the treatment of diffusion, hy- 
drogen, hélium and the metals. But in the most advanced 
version of the code, the diffusion of a dozen species is con- 
sidered: the species involved in the nuclear reaction network 
(except Li, Be and B), Fe and a fictitious species gathering 
ail the other ones. The mean degree of ionization of each 
species is computed and the mean charges of the ions are 
used in the computation of the diffusion coefficients. 

Parametric turbul ent diffusion has also been implemented 
(Miglio etalj|2007d) and in on e version of the code, radia- 
tive forces from the OP project dSeatonli2005) are availa ble 
dBourge and AlecianLl2006al : lBourge et al.ll2006bLl2007h . 
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3.7 Mass loss 



4.2 Chemical évolution 



The inclusion of a diffusion process without mass loss leads 
to overabundances of hydrogen in the external layers. For 
solar models or stellar models with convective envelopes, 
this effect is hidden by the mixture due to convection. But for 
more massive stars, specially when radiative forces are taken 
into account, it turns out that a slight mass loss, compatible 
with admitted values, (of the order of 10~^ M0/yr for a 10 
Mq star) is necessary to prevent the building of unrealistic 
abundances in the external layers of the star and has been 
implemented (Bourge et al., 2007). Work is in progress to 
cope with higher mass losses encountered in more massive 
stars. 



4 Discretization 

The physical quantities describing the star are recorded at 
a number of discrète points. Though their numbers can vary 
during the course of évolution (by addition or deletion), thèse 
points are linked to a given material élément of the star, in 
other words, as long as point Pj^ is not deleted, the value 
of niji is preserved. In fact, the fundamental quantities used 
to characterize our lagrangian grid are not the otj. but the 
masses of the shells between the grid points — m^,, for 
the sake of précision in the external layers. 



4.1 Structure 

At each epoch t, a model is computed by solving the follow- 
ing différence équations. 



1 



dP\ /'dP\ 



m+i-m = x(Pyt + Pi+i)(V4+i -Vi), 



(46) 



(47) 



L^+i -L^ = -(e^ + e^i + e^+i + eg,i+i)('«i+i -'«i), (48) 



71 = 



dT_ 
1^ 



dT_ 
dr 



k+\. 



{rk+\~rt), 



(49) 



with V = 47rP /3. In the energy équation (1481) , the e,, terms 
involve différences between quantities at epoch t and t — 5t. 
Omitting index k. 
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3êT[2{X -X'') +0-75 (y - y") 
+0.5(Z-Z*')]}/5f, (50) 
where upper index refers to the model at epoch t — 5t. 



We update the chemical composition in a two-step process. 
The change of composition due to diffusion is first com- 
puted. In a second step, the change of composition due to 
nuclear burning and mixing are computed simultaneously. 
This two-step process reflects the history of the development 
of the code (diffusion was not included in the first versions) 
and will be corrected in future versions. Fortunately, the dif- 
fusion process turns out to be fast only in the external layers 
of the star where there is no nuclear burning. The code is not 
able to follow the évolution of Li, Be and B. 



Diffusion. In certain circumstances, the diffusion is fast 
enough in the external layers, so that the time scale of the 
diffusion is smaller than any reasonable timestep (stiff dif- 
ferential équations). In such cases, numerical stability of the 
discrète scheme requires that the second members of the 
évolution équations Q be calculated using the updated val- 
ues. 

With index i referring to the chemical species, index k to 
the grid point and upper index to the previous epoch t ~ 5t, 
our discrète équations may be written (omitting unimportant 
détails linked to the non equidistance of the mesh points) 

^^^^ = -{{4nr^pXm)k+i/2 - (47rr2pX,w,),_i/2} 

/(»ÎA:+l/2-»î/t-l/2), (51) 
Wa+l/2 = {j^'^hj,k+l/2i^j.k+l ~^jk) 

i 

+a/',;i+i/2(ln/'A:+i -InPi) 

+ar,*:+i/2(ln71+i -ln7*;)|/(''A-+i - rj,), (52) 

where the fractional values of the spatial index have obvious 
meaning: 



rk+\l2 = {rk + rk+\)/2, ... 



(53) 



Nuclear reactions. The life-time of the différent nuclear spe- 
cies are of very différent orders of magnitude, with some 
of them very much shorter than the timestep. We apply the 
same technique as described in the previous paragraph to 
this stiff differential System. As the life-time of hydrogen is 
of roughly the same order as the duration of the main sé- 
quence, it is possible (and that gives better approximation) 
to use an average abundance of hydrogen X = {X +X^)/2 
in the second member. This possibility has not been used in 
Clés. 



5t 



ijk ) 



(54) 



where the indices have the same meaning as above, except 
that index refers now to the abundances obtained after the 
diffusion step. 
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In a convective zone, the Xn^ are equal for ail points k 
in the zone and équation ( 1541) must be replaced by the dis- 
cretized form of (Hll. However a simple and direct discretiza- 
tion would produce unphysical discontinuities when a con- 
vective zone recèdes by more than one spatial interval in 
one timestep. A dedicated algorithm has been implemented 
to produce a continuons chemical profile in this case. 



5 Implementation 

5.1 The gênerai structure 

IStancliffj ( l2006l) distinguishes three différent stratégies for a 
stellar évolution code to advance by one timestep. Our code 
falls in his partially simultaneous category. A timestep be- 
gins with an updating of the chemical composition followed 
by solving the structure équations. Thèse two processes are 
then repeated until convergence. 

Updating the chemical composition. The discretized équa- 
tions of évolution are nonlinear. They are solved with the 
Newton-Raphson algorithm. We use the composition of the 
previous model as a first guess. 

Solving the structure équations. The structure équations are 
solved by the Newton-Raphson algorithm. When it gives 
signs of difficult convergence, we apply only a fraction of the 
correction suggested by the algorithm and when this recipe 
fails, we adopt a smaller timestep. In the présent state of the 
code, we use the last computed model as a first guess to start 
the itération. We could probably improve the rate of conver- 
gence by using an extrapolation of the last two computed 
models as first guess. 



5.2 Interpolation in tables 

To guarantee the continuity of interpolated functions as well 
as their first derivatives, we use cubic splines. In 1-D, con- 
tinuity of the second derivative can be achieved with cubic 
splines. However, this results in non-local and sometimes 
unphysical behaviour of the interpolating function. We pre- 
fer to be able to impose the values of the first derivative at 
the grid points, values that we compute locally (involving 
only three points), and therefore we sacrifice the continuity 
of the second derivative. This interpolation strategy can be 
extended to more than one dimension in several ways. When 
interpolating physical data, affected by their inhérent uncer- 
tainties, it seems that the choice between the différent vari- 
ants is unimportant. The interpolation is faster if the deriva- 
tives are stored with the values, but the tables occupy more 
space in memory. Both stratégies have been used in Clés, we 
tend now to favour the gain in Computing time. 



5.3 Choice of the grid 

The grid of points of the model is adapted during the évo- 
lution in such a way as to limit the variations of différent 
physical quantifies from one point to ifs neighbours. The 
controlled variables are r, m, P and T . The grid of one model 
is obtained from the grid of the previous one, with addition 
or deletion of points so that the following rules are satis- 
fied: Ar < S.IQ-^R, Am < 5.10"^M, AP < S.IQ-^P and 
AT < 10~^. A simple command allows the user to impose 
stricter or looser rules for the grid choice. 

With the default rules, a 2 model typically starts on 
the Hayashi track with 700 points, reaches the zero-age main 
séquence with 1 150 points and keeps about the same number 
of points along the main séquence with local additions and 
deletions. 



5.4 Choice of the timestep 

The timestep is chosen small enough to ensure that selected 
physical quantifies do not vary too much from one model to 
the next one. The controlled variables are the local values 
of r, P and L and the central hydrogen abundance Xc. The 
allowed maximum variations dépend on the évolution stage 
and are read from a table. This strategy gives the flexibility 
to adopt less severe criteria for the pre-main séquence com- 
putation for instance. However, when the Newton-Raphson 
algorithm fails to converge, the timestep is halved. The user 
has the possibility to influence, in a certain measure, the 
choice of the timestep. 

With no intervention of the user, a 2 Mq model évolu- 
tion needs 125 steps for the pre-main séquence, 75 steps for 
the main séquence and 20 steps more for the second gravita- 
tional contraction. 

However, for évolution séquences computed with diffu- 
sion, the timestep is effectively controlled by the necessity 
of convergence of the diffusion algorithm. 



5.5 First model 

When the luminosity of the star is mainly powered by the 
nuclear reactions, it is enough to give the chemical composi- 
tion as a function of the mass, X, (m) to characterize a model. 
The good model will not be far from the model obtained by 
solving the structure équation with = 0. This is not true 
for our first model on the Hayashi séquence. To characterize 
it completely, we have to impose, for instance, ifs luminosity 
and assume that it is contracting in a homologous way. 



6 The use of Clés 

There are différent levels of sophistication in the use of Clés. 
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When it starts, Clés reads a command file, i.e. a succes- 
sion of commands composed of a keyword followed by pa- 
rameters, instructing the program of the task it has to per- 
form: mass of the star, chemical composition, EOS table, 
opacity table, convection parameter,. . . Clés supplies default 
values for the parameters not specified by the user. The basic 
user must just know how to write a command file to express 
his requirements. This is well explained in a User 's Guide 
of ten or so pages. 

The next step in the use of Clés consists in providing the 
program with new tables, such as opacities for a new métal 
mixture. Tools are available to help the user in this task. 

The more informed user may need to test conditions dur- 
ing the évolution and interact with the program. A subrou- 
tine called clesuser, which does nothing in the standard 
version of the program, can be customized to meet the need 
of the user. It is called by the program before each timestep 
and gives the user the opportunity to read or modify the last 
computed model or the parameters used for the computation. 

At last, a few users have implemented in their copy of the 
code the features needed for their particular research. This is 
facilitated by the modular structure of the code. 



Table 1 Parameters of calibrated solar models. The normalized radius 
at the base of the convective envelope icz and the surface hehum abun- 
dance (Ys) are aiso reported. The values of thèse par ameters obtained 
by helioseismology are respectively 0.7 13 ± 0.001 dBa s u and AntiaL 
,1997.) and 0.245 ±0.005 ( Boothroyd and Sackmamll2003D . 



Model 


a Xo 


Zo 


Tcz 


Y, 


OPAL 


1.805 0.7043 


0.0201 


0.7156 


0.2450 


OP 


1.818 0.7056 


0.0201 


0.7133 


0.2442 



OP 
OPAL 



0.005 



7 Calibration of solar models 

The solar models are calibrated by adjusting the mixing- 
length parameter a, the initial mass fraction Xq of hydrogen 
and the initial mass fraction Zq of heav y éléments in qrder t o 
reproduce, at the solar âge (4.57 10^ yr. lBahcall et al."l995l): 
the o bserved luminosity (3.8 42 10^^ erg s~\ Bahcall et_al, 
2001 , radius (6.9599 lO'" cm,'AlWl973') and ratio of heavy 



r/R 

Fig. 1 Relative différences in the squared sound speed between the 
calibrated models and the Sun ( Basu et al.. 2000.) . 



relative différences in t he squared sound speed inferred from 
helioseismic inversion (iBasu et al .1120001) is shown in Fig.[T] 



8 Discussion 



éléme nts to hydrogen in the photosphère (Z/X=0.0245. lGrevesse and Noelsl 

Il993h . In the solar calibration we adopt a spatial and tempo- 

ral discretization finer than the standard one (see Sec.|5]l: the 
number of mesh-points in a solar model is typically ~ 2500 
and ~ 350 timesteps are needed to reach the solar âge (in- 
cluding the pre-main séquence phase). 

As the values of the calibrated parameters dépend on the 
choice of the physical inputs, as an example we here présent 
two solar models co mputed with différent opacity tables: 
OPAL dîglesias and R ogersl 11 9961) and the récent OP opac- 
ities (Seaton, 2005) both complemented at low température 
with JFfirguson_et_ai. (2005) opacities. The remaining com- 



mon physical in puts are: the mixing - length theory for con- 
vection, GN93 (Grevesse and Noelsl 119931) métal mixture, 
NACRE+Formicola et al. (2004) nuclear reaction rates and 



From the standard version 18 of Clés, a few users/developers 
working in close collaboration in the Asteroseismology Group 
of Liège have implemented a number of désirable features in 
the code. At the moment, there is no single copy of the code 
implementing ail the features described above. Our first goal 
is now to gather the most significant developments already 
achieved and a few more in a unique version 19. And we 
hope that this cycle of development will be repeated to cre- 
ate new versions. 

This strategy of development was made possible because 
the developers have been working in close collaboration on 



the OPAL 2001 équation of state dRogers and Nayfonovil2002h .thp same site. In less favourable circumstances, it would re- 
The outer boundary condition is give n by a radiative gray quire a good documentation. Unfortunately, this is not the 
Eddington atmosphère computed using iFerguson et al.l (12005) case, except for a user's manua l. We are still far from the 
opacities. The models include microscopic diffusion of H, idéal situation described by|^ (|2006|) but we will make the 
He and Z as in the standard version of the code (see Sec. writing of a good documentation our second priority. 
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The parameters resulting from the calibration are reported 
in Table [T] along with the normalised radius of the convec- 
tive envelope and the hélium abundance in the envelope. The 
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